International Journal of Advances in Engineering & Technology, Jan 2012. 

©IJAET ISSN: 2231-1963 



Analysis and Simulation of Series Facts Devices to 
Minimize Transmission Loss and Generation Cost 

M. Balasubba Reddy , Y. P. Obulesh and S. Sivanaga RajiT 

! Deptt. of Electrical Engg., Prakasam Engineering College, Kandukur, A. P., India. 

2 Deptt. of Electrical Engg., L.B.R. College of Engineering, Mylavaram, A.P., India. 

Deptt. of Electrical Engg., J.N.T.U. College of Engineering, Kakinada, A. P., India. 



Abstract 

This paper presents the severity of over load index (SOL) technique to find the optimal location of FACTS 
controller to achieve optimal power flow solution. The proposed method estimates the SOL index of 
transmission line to locate the FACTS device. Main objective of optimal power flow is to minimize the overall 
cost function which includes total active and reactive power production cost under constraints. Among various 
FACTS devices thyristor controlled phase shifters (TCPST), thyristor controlled series compensator (TCSC) 
and static synchronous series compensators (SSSC) are considered to control active power flow through 
transmission line and also to reduce active power loss. Different operating conditions of the power system are 
considered for finding the optimal location of FACTS controllers. The proposed technique is an effective method 
for the optimal location of FACTS controllers. 
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I. Introduction 

In today's highly complex and interconnected power systems, there is a great need to improve electric 
power utilization while still maintaining reliability and security. Demand of electric power is 
continuously rising at a very high rate due to rapid industrial development. To meet this demand it is 
essential to raise the transmitted power along with the existed transmission facilities. The need for the 
power flow control in electrical power system is thus evident. Power flow is a function of 
transmission line impedance, the magnitude of the sending end and receiving end voltages and the 
phase angle between voltages. By controlling one or a combination of the power flow arrangements, it 
is possible to control the active, as well as, the reactive power flow in the transmission line [1]. The 
FACTS devices are capable of changing the system parameters in a fast and effective way. It is 
known that the benefits brought by FACTS devices include improvement of system stability, 
enhancement of system reliability, and reduction of operation and transmission investment cost [2]. 
In previous work, researches concentrated on locating and sizing of different types of FACTS devices 
in order to maximize the power transfer considering networks with variable loads. The problem was 
formulated as an optimization problem and was solved using different methods such as using iterative 
techniques, MATLAB optimization routines or Genetic Algorithm (GA) [3-5]. Optimization problem 
with the objective of minimizing the generating cost in a network with unchanged loads. The problem 
is solved using OPF algorithm using NR method. Which is fast and simple compared with 
conventional technique and also give promising results. Improvements of results with FACTS 
devices is compared with convention N-R OPF method without FACTS devices. OPF is a very large, 
non-linear mathematical programming problem, the main purpose of OPF is to determine the optimal 
operation state of a power system while meeting some specified constraints. Since the OPF solution 
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was introduced by squires [6]. The focus in this paper lies on thyristor controlled phase shifting 
transformers (TCPST), thyristor controlled series compensators (TCSC) and static synchronous series 
compensators (SSSC). These facts devices TCPST, TCSC and SSSC by controlling the power flows 
in the network can help to reduce the flows in heavily loaded lines, to reduce the system loss, reduced 
cost of production and fulfilled contractual requirements [7]. Location of FACTS devices in the power 
system are obtained on the basis of static and/or dynamic performance. There are several methods for 
finding the optimal location of FACTS devices but this paper introduces severity of overload index 
technique gives optimal location by estimating overload index of each transmission line. The 
proposed algorithm is tested using the IEEE 5 bus and IEEE 30 bus system and the results are 
presented. 

The description of remaining sections of this paper is as follows. Section 2 introduces OPF without 
FACTS devices. Modelling of FACTS devices is described in section 3. In section 4 problem 
formulations is presented. The experimental results on the IEEE5 bus and IEEE30 bus systems are 
presented in section 5. Finally section 6 summarizes the main conclusion and future scope of the 
paper. 

II. OPF without Facts Devices 

The objective of active power optimization is to minimize production cost while observing the 
transmission line and generation active and reactive power limits. 
Minimize ™ 

F T = Zj C i\ P Gi) (1) 



;=i 



Subjected to 



IX-2X-^= 



(2) 

i=\ k=\ 

p < p max 

r L ~ ^ (3) 

pmin < p < pmax 
r Gi — r Gi — r Gi (4) 

Where n is the number of system buses and number of generating units respectively. C, [P Gi ) is 

production cost of the unit at i' bus, F T is the total production cost of m generators, P^ m & P^ dX are 

minimum and maximum active power limits of the unit at i bus. P Dk is the active power load at bus 

k, P L is the network active power loss, P b Pj max are the active power flow and its limit on line 1. 
The augmented lagrangian is, [18] 

(n m ' 

,=i ; (5) 

f>>, -pr)+t\Mr(p c T-p Gl )+Mr(p Gl -p G r x )} 

A is for power balance equation. 

//™ n and// ; max are lower and upper active power limits of unit at i bus. 

jU[ is for active power flow limit on line 1. 

N t is the number of transmission line flow violations 

III. Modelling of Facts Devices 

In this paper three facts devices are used for comparing the results these are thyristor controlled phase 
shifting transformers (TCPST), thyristor controlled series compensators (TCSC) and static 
synchronous series compensators (SSSC). Modelling of these devices is given in the following 
sections. [20] 
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3.1 Modelling of TCPST 

TCPST can be modelled by a phase shifting transformer with control parameter (j). The power flow 
equations of the line can be derived as follows [15]& [17] 
P k =V k 2 G- V k V m [G Cos{6 k - 9 m - </>) + B Sin{d k - 8 m - </>)] 

Q k =-V k 2 B-V k V m [GSir{e k -d m -$)-BCo{d k -6 m -<?>)] 



P m =y,nG-V„y k [GCo{O m -9 k -0)+BSir(0 m -e k -+)] 

Q m =-v m 2 B-v m v k [GSir(e k -e m -<p)+BCo{e k -e m -*)] 

The injected real and reactive power of TCPST at bus and bus are as follows. 
P ks = -V t 2 G 9 tan^-nV, tanfc SinS kl -B kl CosS kl ] 
P ls = -V k V, tm4G kl SinS kI +B kl CosS kl ] 



(6) 

(7) 
(8) 
(9) 

(10) 
(11) 



3.2 Modelling of TCSC 

It was derived by examining the voltages and currents in the TCSC circuit under the full range of 
operating conditions. The basic equation is, 



Z TCSC =R TCSC +jX 7 



v 

_ ' TCSC 



l TCSC 



(12) 



Where, V TCSC is the fundamental frequency voltage across the TCSC module, 



"TCSC 



is the fundamental frequency line current and Z TCSC is the TCSC impedance. 



Effective TCSC reactance X TCSC with respect to alpha is, 
x tcsc ( a ) = ~ x c + c \ ( 2 (^ -a)+ Sin{2{x - a))) - 
C 2 Cos 2 {x - a){co(ta.n{(o(ff - a))) - tan(^ - a)) 



(13) 



3.3 Modelling of SSSC 

The SSSC OPF model presented in this section enables very flexible and reliable power system 
optimization studies to be carried out. The flexibility stems from the generality of the SSSC model 
and the robustness from the strong convergence exhibited by the OPF solution using Newton's 
method. The SSSC model may be set to control active powers as well as nodal voltage magnitude, at 
either the sending or the receiving end bus. The mathematical modelling will be obtained from [8] and 
that was taken as a reference and implemented in the OPF by using the Lagrangian multiplier concept. 
A SSSC [9] [10] usually consists of a coupling transformer, an inverter and a capacitor. As shown in 
Fig. 1 the SSSC is series connected with a transmission line through the coupling transformer. 

Bus i Vsg Busj 



Ii 



Vi 



TJHMLT 

room 



Delink 




Ij 



Fig. 1 SSSC basic operating model 



It is assumed here that the transmission line is series connected via the SSSC bus j. The active and 
reactive power flows of the SSSC branch i-j entering the bus / are equal to the sending end active and 
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reactive power flows of the transmission line, respectively. In principle, the SSSC can generate and 
insert a series voltage, which can be regulated to change the impedance (more precisely reactance) of 
the transmission line. In this way, the power flow of the transmission line or the voltage of the bus, 
which the SSSC is connected with, can be controlled. 

An equivalent circuit of the SSSC as shown in Fig. 1 can be derived based on the operation principle 
of the SSSC. In the equivalent, the SSSC is represented by a voltage source V^in series with 

transformer impedance. In the practical operation of the SSSC, V Vf can be regulated to control the 

power flow of line i — j or the voltage at bus /or j [11]. 

In the equivalent circuit, V se =V se Z.B se ,V i =V i Z.B i ,V se =V j Z.0jthea the bus power flow 

constraints of the SSSC are, 

P,j =V-Sn -V^g ycosffy +b v sm9y)-Vy e (gyCos{9.-9 se )+by sin(0. -&J (14) 

Qtj =-v-h -vy^y sin^ -by cosOyj-vy^gy Mo.-eJ-b^ cosfe. -ejj (is) 
p ]l =vjg ]r yy,(g 1] oo^ jl +b l] ^e j )+vy s j^ l] co^ ] -e s },+b l] ^e ] -e s ]) (16) 

Qji =- v >j J -Wjigy sin^ -by cos0..)+V.vJg, sin(^. -O se )-by cos(^ -9 J (17) 
where, 

gij + jby =^,g H = gij A =b ipgjj = gl] , bjj =b tj The operating constraint of the SSSC (active 
power exchange via the DC link) is, 

PE = Rc(vj;) = (18) 

where, 
Re (Yselji) = ViVseidij cos(0i - se ) - b tj sin(0; - d se ))VjV se {g t j sm(6j - Q se ) - b tj cos(fy - 

e se j) (19) 

EiiidjQii Eii+iOii 

Bus I Zse Vse Busj 




efVse Iiil = 



IV. 



Fig. 2 SSSC Equivalent circuit 

OPF Formulation with SSSC Lagrangian Function 



Based on the above equation the Lagrangian function for the SSSC may be written as, 

L{x,A) = f(p g )+A'h(p g ,V,9,S cr ,V cr ) (20) 

In this expression, f\P j is the objective function to be optimized /z(P ,V,0,V CR ,S CR j; represents 
the power flow equations; x is the vector of state variables; A is the vector of Lagrange multipliers for 
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equality constraints; and P ,V and Owe, the active power generation, voltage magnitude, and 
voltage phase angle, respectively. The SSSC control variables sxed CR andy cR . The inequality 
constraints g\P 'V,0,V CR ,S CR )< are not shown in Eqn. (3.81) because it is added only to l(x,A) 

when there are variables outside limits. [19] 

The Lagrangian function, L km (x,X) corresponding to the power flow mismatch equations at buses k 

and m, is given by the following equation, 

Lkm( x >ty = A-pkVk + Pdk ~ Pgk) + ^qk\Qk + Qdk ~ Qgk) + Vi'™ + ^dm ~ Pgm) + 

AqmiQm + Qdm ~ Qgm) ( 2l ) 

In this expression P dk ,P d ,Q„ and Q dm are the active and reactive power loads at buses k and m; 

Pvk'Pgm->Qgk an d Qem are tne scheduled active power generations at buses k and m; and 
A k ,A ,A . and A are Lagrange multipliers at buses k and m. The vector of state variables x is 
[V S] < where V and ^ include both nodal voltages and SSSC voltage sources. 

4.1 SSSC Power Flow Constraints 

The power injected at bus m by the SSSC can be formulated as a power flow constraint in the branch 
connecting buses m and / . We may write 

'-'ml \ X ->^) = Apml Vml ~ "specified ) + \ ml \*£ml ~ *Z specified I (22) 

Where, A Pml is the Lagrange multiplier associated with the active power injection at bus m; and 

Specified ^ s tne specified active power leaving bus m. 

In conventional OPF formulations, such constraints are enforced only if power flow limits have been 
exceeded. However, in this particular application this constraint may remain active throughout the 
iterative solution [12]. The SSSC Lagrangian function comprising the individual contributions 
presented above is as follows 

L sssc {x,A) = L km {x,A)+ L se {x, A) + L ml {x,A) (23) 

4.2 Optimal setting of SSSC Parameters 

Basically, the SSSC has one voltage source inverter (VSI) having a dc storage capacitor. It is 
connected to the system through a coupling transformer. In this study, the series compensation 
^ facts =A£7 SSC is employed [12]. 
The injected currents at bus /and bus j can be expressed as follows, 

At/ 



A/. . 

i, send 



Z tj (24) 



A/. d =-^ 

Z, (25) 

The cost function for SSSC can be expressed as follows, 

C Jsse = 0.00 15s 2 -0.5130^ + 133. 15 $/hr (26) 

where, C is the cost function of the SSSC in $/hr and V is the operating range of the FACTS 

devices in MVAr. 

V. Severity of over Load ability index (SOL) Computation 

The location of the FACTS devices in this work is decided based on the severity of the overloading of 

that particular branch in which the device is incorporated. The process of ranking the branches based 

on their load ability in the order of their severity involves the following steps. 

Stepl : Establish the criterion to be considered in formulating the ranking 

Step2: For the criterion established in (Step 1), define a scalar mathematical function which has a 
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large value of branch load that which stress the system relative to that criterion, and a small value for 
those which do not; this function is called a "SOL index." 

The SOL index is such that contingencies resulting in system conditions yielding large valued over 
load indices are considered more severe than system conditions with smaller over load indices [16]. In 
the overload ranker, the SOL index is defined as, 

2 



SOL 



n 

i 

i = \ 



( 



Pi 

P i , m a x J 



(27) 



where, 

P is the real power flow in line "/", 

P t max is the maximum of active power transfer over the i' line and 
' n ' is the set of monitored lines contributing to SOL. 
5.1 Calculation of SOL for IEEE 5 Bus system 

Table 1: SOL index of all buses by running the general OPF for IEEE 5 bus system 



Bus No., Node 
No. 


SOL index of 
each bus 


Ranking 


[3] 


0.5S12 


1 


[4] 


0.5310 


2 


[5] 


0.3285 


3 



As compared the above SOL-indices for the IEEE 5 bus system among the 3 load buses (3, 4, 5) the 
bus 3 is having the maximum SOL index, it is considered to be the critical bus. Hence line indices will 
provide accurate information with regard to the stability condition of the lines. 
5.2 Calculation of SOL for IEEE 30 Bus system 

Table 2: SOL-indices by running the general OPF of maximum loaded buses in IEEE 30 bus system 



[Branch 

Number 


SOL indices 

of different 
branches 


Ranking 


[30] 


7776 


<> 


[24] 


05672 


4 


[29] 


0.SS73 


1 


[28] 


0.7486 


3 


[26] 


0.5491 


5 



As we considered the SOL-index table of the IEEE 30 bus system there will be the 5 load buses (24, 
26, 28, 29, 30) with the bus (29) is having the maximum load ability, it is considered to be the critical 
bus. The branch connected to that particular weakest or critical bus will be the optimal location for the 
FACTS device to be placed. Hence the branch [13]-[14] is chosen to be the optimal location in the 
IEEE 30 bus case. 

VI. Simulation Results 

The proposed methodology of active power optimal power flow of series FACTS devices as TCPST, 
TCSC and SSSC for transmission network is implemented using MATLAB on the two test systems 
viz., IEEE 5-bus and IEEE-30-bus test systems. The cost constants and their typical values for this 
problem are energy cost (a=60 Rs/hr, b=200 Rs/MW/hr, c=140 Rs/ ^1^2 /hr). The payback period 

assumed is five years with depreciation factor (qA as 0.1. 

6.1 Optimal power flow for IEEE 5 bus system without and with FACTS devices 

The 5-bus test system used to illustrate the use of the conventional power flow Newton- Raphson 
method and is also used to illustrate the use of the OPF and associated data. The maximum and 
minimum voltage magnitude limits at all buses are taken to be 0.9p.u. and 1.1 p.u., respectively, 
except at bus 1 (Slack bus), where the maximum limit is set at 1.5 p.u. 
The cost coefficients of the two generating buses are taken to be 

a = 2700 RslMw 2 Ihr 

b = 153 Rs I 'Mw I hr and, 

c = 0ASrs/hr. 



468 



Vol. 2, Issue 1, pp. 463-473 



International Journal of Advances in Engineering & Technology, Jan 2012. 

©IJAET ISSN: 2231-1963 

The maximum and minimum generator active power limits are set at 200Mw andlOMw, respectively, 
whereas the maximum and minimum reactive power limits are set at 300MVAR and -300MVAR, 
respectively. 

Table 3 Comparison of OPF solution without and with FACTS devices for IEEE 5 bus system 



Without device 


TCPST 


TCSC 


SSSC 


VW-u. 


VbusP.U. 


V bu3 p.u. 


V bsu p.u. 


1.060 


1.06 


1.060 


1.060 


1.100 


1.012 


1.035 


1.049 


1.078 


1.028 


1.042 


1.048 


1.077 


1.027 


1.028 


1.036 


1.072 


1.020 


1.032 


1.049 




1.009 


1.018 


1.027 



Table 4 Comparison of Active power and Active power generation cost of different cases for IEEE 5 bus 

system 



Case type 


Piwi (MW) 


Petmstm (Rs/hr) 


Without device 


3.55 


34,046 


TCPST 


3.45 


33,874 


TCSC 


3.41 


33,558 


SSSC 


3.21 


33,366 



I 



— ■-W/o device 
-♦■-TCPST 
—•—TCSC 
-*-SSSC 




3 3.5 

No. of Buses 



4.5 5 5.5 6 



Fig.3 Comparison of without devices and with TCPST, TCSC and SSSC voltages magnitudes at each bus 
in IEEE 5 Bus system 



© - 5 



<5-. 




Tv o of test cases 

Fig.4 Comparison of Active power loss and % of P , without and with FACTS devices 



469 



Vol. 2, Issue 1, pp. 463-473 



International Journal of Advances in Engineering & Technology, Jan 2012. 

©IJAET ISSN: 2231-1963 

6.2 Optimal power flow for IEEE 30 bus system without and with FACTS devices 

The IEEE 30-bus test system is used to study the impact of the FACTS devices on the network. The 
optimal location will be analyzed by using the SOL index from the Table 2 and then FACTS devices 
are added in series with transmission line (28-29), and the dummy bus (31) is added to enable such a 
connection to take place. 

As compared the voltage magnitudes without and with FACTS devices the voltage profile of the 
SSSC is improved slightly when compared with TCPST and TCSC, for both the IEEE Test systems. 
By comparing the SOL-index under normal situation the optimal location of the FACTS device is 
decided. Hence for the IEEE 30 Bus system (28-29) is the optimally decided branches for the FACTS 
devices to be incorporated in the electrical power system. 



Table 5 The active power and reactive power of different buses for IEEE 30 bus system 


Case type 


Active power loss (MW) 


Reactive power loss (MVAr) 


Without Device 


18.58 


52.73 


With TCPST 


18.49 


46.41 


With TCSC 


18.37 


49.23 


With SSSC 


18.02 


42.45 



Table 6 The initial and final costs of active power generation at different buses 



Case type 



Without Device 



With TCPST 



With TCSC 



With SSSC 



IEEE 5 bus system 



Gbh, initial cd& 
(Rs/hr) 



35,000 



35,000 



35,000 



35,000 



Gbh, final ct>st 
(Rs/hr) 



34,046 



33,874 



33,558 



33,366 



IEEE 30 bus system 



Gbh, initial cds; 
(Rs/hr) 



36,900 



36,900 



36,900 



36,900 



Gbh.. final tost 
(Rs/hr) 



36,765 



35,325 



33,795 



33,075 



Furthermore, with SSSC the generation cost is reduced to 680 Rs /hr in 5 bus system and 
3825 Rs/hrin 30 bus system respectively when compared to TCPST , TCSC and with the base case i.e. 
4.66% reduction in the active power generation cost compared to (2.76% , 3.25% and 4.12% for the 
without FACTS .with TCPST and TCSC respectively). 

Table 7 Comparison of Total Active power generation costs in both the test cases 



Case Type 
device 


IEEE 5 bus system 


IEEE 30 bus system 


Jmii jira ( MW ) 


'WwsWfct-a; (Rs/hr) 


(MW) 


F FfTf*MreW«t-5 : . (Rs/hr) 


Without device 


3.55 


34,046 


18.58 


36,765 


TCPST 


3.45 


33,874 


18.49 


35,325 


TCSC 


3.41 


33,558 


18.37 


33,795 


SSSC 


3.21 


33,366 


18.02 


33,075 




No. of Buses 



Fig. 5 Comparison of without and with TCPST, TCSC and SSSC voltages at different buses in IEEE 30 Bus 

system 
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Fig. 6 Comparison of Active power generation cost (Pgi, cost) in case of without device & with 
TCPST, TCSC and SSSC at different buses in IEEE 30 Bus system 

Table 8 Comparison of the costs of different FACTS devices in IEEE 5 bus & IEEE 30 bus system 



Type of Device 


Cost of the FACTS devices [(Rs/hr) for TCSC & SSSC & (Rs) for TCPST] 


IEEE 5 bus system 


IEEE 30 bus system 


TCPST 


1,86,075 Rs 


3,66,075 Rs 


TCSC 


6,938.99 Rs/hr 


17,600 Rs/hr 


SSSC 


5,952.47 Rs/hr 


13,452 Rs/hr 



From Tables (3) and (6) it was clear that the SSSC is improving the overall system performance in 
both the IEEE test systems when compared with the TCPST, TCSC and base cases. Being having the 
separate DC link for the SSSC it is capable of generating the reactive power without any external 
source and hence capable of controlling the reactive power flow of the system like UPFC, instead of 
controlling only active power as TCSC. The Table 5.18 will gives the costs of the FACTS devices 
among which the SSSC is possessing the less cost comparable to the other devices. Hence it was 
concluded that the SSSC is more suitable for the power flow control and optimizing techniques. 

VII. Conclusions 

Simulation studies using MATLAB programming code, on IEEE 5 and IEEE 30 bus system are 
presented to illustrate the methodology and to demonstrate the benefit of the proposed approachln this 
paper the optimal power flow analysis with the inclusion of TCPST, TCSC and SSSC has been done. 
Newton-Raphson method used in polar co-ordinate form is effectively applied to solve the Optimal 
power flow equations of IEEE 5-Bus & IEEE 30-Bus systems which differ from each other in size & 
degree of operational complexity. 

Themodeling of Series FACTS devices are incorporated into an existing Newton-Raphson load flow 
algorithm, which is capable of solving large power networks very reliably. It shows that Series 
FACTS environment can be set to control active power, active power loss and voltage magnitude 
simultaneously. The Severity of Overloading index is used to find location of the series FACTS 
devices that they are used in the network. 

The proposed algorithms were implemented to find out the proper setting and installation cost of the 
TCPST, TCSC and SSSC in IEEE-5 bus & IEEE-30 bus test systems. By comparing the results, it is 
observed that SSSC is more effective than TCSC and TCPST in terms of voltage regulation, power 
loss reduction, minimization of active power generation cost and improving the active power flow. 
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